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Recent numerical simulations have found that the Cauchy horizon in- 
side spherical charged black holes, when perturbed nonlinear ly by a self- 
gravitating, minimally-coupled, massless, spherically-symmetric scalar field, 
turns into a null weak singularity which focuses monotonically to r = at 
t ■ late times, where the singularity becomes spacelike. 

■ Our main objective is to study this spacelike singularity. We study ana- 

lytically the spherically-symmetric Einstein-Maxwell-scalar equations asymp- 
totically near the singularity. We obtain a series-expansion solution for the 
metric functions and for the scalar field near r = under the simplifying 
assumption of homogeneity. Namely, we neglect spatial derivatives and keep 
only temporal derivatives. We find that there indeed exists a generic space- 
like singularity solution for these equations (in the sense that the solution 
depends on enough free parameters), with similar properties to those found 
in the numerical simulations. This singularity is strong in the Tipler sense, 
namely, every extended object would inevitably be crushed to zero volume. 
In this sense this is a similar singularity to the spacelike singularity inside un- 
C***~ . charged spherical black holes. On the other hand, there are some important 

differences between the two cases. Our model can also be extended to the 



o 
in 



5-H 



more general inhomogeneous case. 

The question of whether the same kind of singularity evolves in more real- 
istic models (of a spinning black hole coupled to gravitational perturbations) 
is still an open question. 

I. INTRODUCTION 

The singularity theorems of Hawking and Penrose [U predict the occurrence of spacetime 
singularities inside black holes under very plausible assumptions. However, they tell us 
nothing about the geometrical and physical nature and properties of these singularities. 
The Kerr black hole, which describes the exterior of realistic black holes, fails to describe 
adequately their interiors. In the Kerr solution the singularity has the shape of a ring, and it 
is timelike. Consequently, it violates the strong cosmic censorship hypothesis of Penrose [pf . 
However, the Kerr solution is unstable, in the sense that the Cauchy horizon is converted 
into a spacetime singularity. This singularity is null and weak || . Namely, infalling extended 
observers would experience only finite tidal distortions at the singularity. Because of the 
complexities involved with the analysis of spinning black holes, a frequently used toy model 
is the spherical charged black hole. Although realistic black holes are not expected to 
be significantly charged, the spherical charged black holes share many common properties 
and similar causal structure with the more realistic spinning black holes. With spherical 
symmetry, non-trivial dynamics can be modeled by a scalar field, which has a spherically 
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symmetric radiative mode. In recent numerical analyses with this toy model (namely, a 
spherical charged black hole perturbed non-linearly by a self-gravitating scalar field) it has 
been shown that the generators of the Cauchy horizon are focused at the late parts of the 
Cauchy horizon fi]||. Namely, the value of the area coordinate r decreases slowly and 
monotonically along the Cauchy horizon, until the Cauchy horizon is completely focused, 
i.e., the area coordinate shrinks to zero value at a finite value of the affine parameter along 
the Cauchy horizon. Then, the singularity becomes spacelike. (The spacelike singularity 
was first found numerically by Gnedin and Gnedin ||.) It is still an open question whether 
in a realistic spinning black hole perturbed non-linearly by a realistic physical field (such as 
gravitational waves) the generators of the Cauchy horizon would also focus completely, and 
create a strong spacelike singularity. In this paper we shall focus on the spacelike singularities 
inside spherical black holes (both uncharged and charged), and study their properties under 
the simplifying assumptions of spherical symmetry and a minimally-coupled massless scalar 
field. We shall also assume homogeneity and power-law behavior of the metric functions, 
which are justified by recent numerical simulations. 

The spacelike singularity inside a spherically symmetric uncharged black hole perturbed 
by a scalar field was studied by Doroshkevich and Novikov 0, who treated the scalar field 
as a linear perturbation and obtained the general solution for the linear scalar field, but did 
not study the backreaction of the scalar field on the geometry. Doroshkevich and Novikov 
found that the scalar field diverges logarithmically with r at the singularity. Therefore, 
the energy density of the scalar field diverges like an inverse power of r, and therefore, 
acting as a source term for the Einstein-Klein-Gordon equations, the scalar field might be 
expected to change the metric functions. The very divergence of the linear perturbation 
analysis indicates that nonlinear effects are expected to be important. More recently, Krori, 
Goswami and Das Purkaystha |5| considered the same problem as Doroshkevich and Novikov 
0, but failed to obtain a general solution. The regular solution obtained by these authors 
depends on just one arbitrary parameter, while a general solution clearly depends on two 
independent arbitrary parameters (see below). Consequently, the correct conclusions to be 
drawn from a linear analysis are the conclusions of Doroshkevich and Novikov [[/J and not 
the conclusions of Krori, Goswami and Das Purkaystha ||. However, with the presence 
of an electric charge a linear perturbation analysis cannot investigate a generic spacelike 
singularity, as this singularity is created by nonlinear effects, i.e., by the nonlinear focusing 
of the generators of the Cauchy horizon. Under linear perturbations there would not exist 
inside a spherical charged black hole a spacelike singularity at all. Therefore, to study 
the spacelike singularity inside a spherical charged black hole with a scalar field one has 
to consider the fully nonlinear case. Nonlinear analysis of spacetime singularities with a 
scalar field in the cosmological case was done by Belinskii and Khalatnikov ||, who found 
that the scalar field destroys the BKL oscillations and that in that case the singularity is 
monotonic. (Interestingly, Belinskii and Khalatnikov showed that if, in addition to the scalar 
field, there were also a vector field, the singularity would again be oscillatory. The analysis 
by Belinskii and Khalatnikov M is markedly different from the analysis we shall present 
below, as we are interested here in spherical symmetry. It turns out that the Kasner-metric 
based Belinskii-Khalatnikov analysis cannot be reconciled with spherical symmetry.) 

The organization of the paper is as follows. In Section || we describe the physical model 
we employ. In Section we discuss the properties of power-law spacelike singularities. 
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In sections [TV] and [V| we study the cases of uncharged and charged spherical black holes, 
correspondingly. We show that indeed there exists a generic solution of the spherically 
symmetric coupled Einstein-Maxwell-Klein-Gordon equations with a spacelike singularity, 
with properties which one actually finds in numerical simulations. This result strengthens 
our confidence in the picture of black hole interiors as described above. In Section [VT] we 
analyze the strength of the singularity we find. We find that it is strong in the sense of 
Tipler (whereas the null Cauchy horizon singularity is weak). We summarize and give some 
concluding remarks in Section [VII| . 



II. THE PHYSICAL MODEL 

We performed numerical simulations of the collapse of a spherically symmetric self- 
gravitating minimally coupled massless scalar field over a pre-existing Reissner-Nordstrom 
black hole. Our code is based on free evolution and on double-null coordinates and is 



described at length in Ref. |I(J. The code is stable and converges with second order. The 
numerical set-up is described in Figure [TJ. Prior to the initial hypersurface the geometry 
is Reissner-Nordstrom. Then, at some advanced time vq the spacetime is perturbed by a 
high-amplitude self-gravitating spherical massless scalar field of squared-sine shape on the 
outgoing section of the initial hypersurface. The null coordinates we use in the numerical 
simulations are linear with the area coordinate r on the initial hypersurface. 

This high amplitude scalar field increases the external mass of the black hole at late 
times by over 10%. We then probe the metric functions and the scalar field in the black 
hole's interior along outgoing (and ingoing) null rays, approaching the spacelike singularity, 
which is located at some finite value of the ingoing null coordinate v*(u p ). Here, u p is the 
value of the outgoing null ray u at which our outgoing null ray is located. Figure [I] shows 
three such rays, denoted by 1,2, and 3. The metric we use for the numerics is 

ds 2 = -F(u,v)dudv + r 2 (u,v)dtt 2 . (1) 

The results we find numerically are as follow: First, the metric function g uv = —F/2 vanishes 
at the singularity, and F decays as a power-law with respect to the area coordinate r as one 
gets closer to the singularity. In addition, we find that the area coordinate r decays like a 
power of — v near the singularity. Numerically, we find in all our simulations this power 
to be very close to |, with a deviation of 1%. Therefore, we would like to show analytically, 
that in a generic solution indeed r oc (if* — v) 1 ^ 2 near the singularity. Our Numerical results 
are described in Figures |2] and [| These results are independent of the precise definition of 
the null coordinates we were using, as any regular gauge transformation which preserves the 
metric-form (|l|) introduces just a multiplicative factor which does not change the power-law 
indices. In addition, we find similar results for all u p 's. That is, the power with which F 
decays to zero at the singularity depends on u p , but only weakly so. The power | we found 
for riv* — v) is retained for all m p 's. 
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Spacelike 




FIG. 1. Spacetime diagram for a spherical charged black hole with a self-gravitating scalar field. 
Singularities are denoted by thick lines. Prior to the initial hypersurface (dashed) the geometry 
is Reissner-Nordstrom. The spacetime is then perturbed by the scalar field. The Cauchy horizon 
is converted into a null singularity. It is still an open question whether there is a continuation of 
the spacetime manifold beyond the null singularity, and if so, what the topology and the geometry 
are. The Cauchy horizon is focused to r = 0, where the singularity is spacelike. Our numerical 
setup does not allow us to investigate the late parts of the spacelike singularity, which are denoted 
by a thick dashed line, as they are located beyond the domain of influence of the characteristic 
hypersurface. We investigate the fields along outgoing null curves, denoted by 1, 2, and 3. 
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FIG. 2. Metric function F = —2g uv as a function of r near the spacelike singularity along an 
outgoing null ray. 

Therefore, we make the simplifying assumption of quasi homogeneity. Namely, we assume 
that only derivatives of the metric functions and of the scalar field normal to the singularity 
are non- vanishing, while derivatives tangent to the singularity are exactly zero on it. We 
thus write the line-element as 
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ds 2 = h{r) dt 2 + f(r) dr 2 + r 2 dtt 2 , 



(2) 



where dfl 2 = d9 2 + sin 2 9 dip 2 is the usual metric on the unit two-sphere. As the singularity 
is spacelike, and t is the spacelike coordinate near the singularity and r is timelike, we are 
going to neglect all derivatives with respect to t, and keep only derivatives with respect to 
r. With these assumptions the t-t, r-r and 9-9 components of the Einstein-Maxwell-scalar 
equations in spherical symmetry are, correspondingly, 



h 



r 2 P 



f'r + f-f 



, q 

f r A 



(3) 



hr 2 



{h'r - hf + h) = $' 2 - / 



(4) 



Af 2 h 2 



(2tifh - 2fh 2 + 2rh"hf - rh' 2 f - rh'f'h) = ^ - y$ /2 . 



(5) 



[The <p-<f) component of the field equations gives again Eq. (|5])]. Here, a prime denotes 
differentiation with respect to r. In addition we also have the Klein-Gordon equation for 
the scalar field = 0, whose first integral reads 



d 



--f(r) sin#, 



(6) 



where d is an integration constant and g is the metric determinant. (Note that because the 
metric determinant has a factor sin 2 9, $'(r) does not depend on 9.) 
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FIG. 3. Area coordinate r as a function of v* — v near the spacelike singularity along an outgoing 
null ray. At small values of r there is growing numerical noise. However, even before the noise 
becomes significant we have over two decades in r and four decades in v* — v of asymptotic behavior. 
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We now use Eq. (U) to eliminate the scalar field from Eqs. (0)-©, and find that 

f'r + f 2 -f + J^{d 2 - q 2 h)=0 (7) 



h'r -hf + h + L (d 2 + q 2 h) = (8) 



In what follows, we shall solve Eqs. (0) and (^j), subject to our conclusions above from the 
numerical simulations. We then use Eq. (||) as a consistency check for our solution. Namely, 
we seek a generic solution for which both metric functions h, f vanish at the spacelike 
singularity. This is motivated by the vanishing of g uv in the double-null metric, from which 
the vanishing of both h and / follows, by merit of the relations 

_ _ du dv 

ytt = a — a ~^t~^t9uv — —^9uv 
at at 

' drS\ 2 ( drA 2 _ ( dr*\ 2 



9rr = f = 9r,r, [ ~T ) = ~9U T" = 2 9 



x dr J \ dr J \ dr J 

Here, the 'tortoise' coordinate r* is defined by g TtTt = —gu, and the null coordinates are 
defined by u = r* — t and v = r* + t. 



III. PROPERTIES OF POWER-LAW SPACELIKE SINGULARITIES 

In this section, we shall study several properties of spacelike singularities, which are 
common to any spatially homogeneous metric of the form @ with h(r), f(r) being powers of 
r. (Note, we still do not assume that f(r), h(r) vanish at the singularity. We shall specialize 
to this case below.) First, we shall give a few examples of known spacelike singularities. 
Then, we shall study the general homogeneous power-law metric (0). 



A. Examples for spacelike singularities 

1. The Schwarzschild singularity 
The Schwarzschild solution assumes near the singularity the form 

ds i ^(^-Y 1 dt 2 - ( 4t") dr 2 + r 2 dn\ (10) 

where t is a spacelike coordinate tangent to the singularity and r is timelike and normal to 
the singularity. This metric obviously satisfies the power- law form of (0). 
The Kretschmann scalar along radial t = const curves is 
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where r is the proper time of the infalling object, and is the proper time at which the object 
arrives at the singularity at r = 0. We now denote by x 4 the coordinate tanget to the world 
line of an infalling observer. Because the Schwarzschild spacetime is vacuous, the projection 
of the Ricci tensor on the world line is identically zero, i.e., -R(4)(4) = R^\^\^) = 0, where 
R( a )(/3) is the af3 tetrad component of the Ricci tensor. These results are independent of the 
mass M of the Schwarzschild black hole. 



2. The Friedmann- Roberts on- Walker cosmology 
The metric for the FRW cosmological model is given by 

ds 2 = -dt 2 + R 2 (t) [~^kri dr " + r * dn ") ■ ( 12 ) 

Let us consider here radiation-dominated perfect fluid with energy momentum tensor T^ v = 
PQfiu + {p + p)u fJi u u , being the four- velocity. We take p = p/3. 
Then, in a co-moving frame near the singularity 

3 1 

RaP^R^ 5 = ^ ^ _ r )4 ' ( 13 ) 



and 



3 1 



%)(4)W = 4^- (14) 



These results are independent of the density p (or the pressure p), and are also independent 
of k, i.e., independent of the spatial curvature. 



3. Kasner solution 



The Kasner solution with a spacelike parameter (similar results are obtained also for the 
Kasner solution with a timelike parameter) is given by the metric 

ds 2 = x 2pi dt 2 - dx 2 - x 2p2 dy 2 - x 2ps dz 2 , (15) 

with pi + p2 + ps = 1 = p\ + p\ + p\. (This solution is anisotropic and is not spherically 
symmetric like the other solutions we consider in this paper.) 
For the Kasner solution one finds 

R aPl sR a ^ 5 = 16p?(l - Pi) , \ w (16) 

and 

%)(4)(r) = 0. (17) 

This last result is obvious from the fact that the Kasner solution represents a vacuum 
spacetime. 
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B. General properties 



In all the above examples, the Kretschmann scalar diverges near the spacelike singularity 
like 1/(t* — t) 4 , and the Ricci tensor, if spacetime is non-vacuum, diverges like 1/t 2 . In this 
section we shall show that these results are general for any metric of the form (fj). Writing 
h(r) = Br m and f(r) = —Ar n , with A, B being constants, a straightforward (but lengthy) 
calculation for the Kretschmann scalar yields 

a g s m 4 — 4m 3 — 2m 3 n + 12m 2 + 4m 2 n + m 2 n 2 + 8n 2 + 16 
Ra ^ sR = 4 (2 + n) 4 (r,-r) 4 • (18) 

In Schwarzschild, m — —1 and n = 1. Substitution recovers the known result for 
Schwarzschild (fLip . This result fll8l) is independent of A,B. This is the analogue of the 
independence of the Schwarzschild Kretschmann scalar of the mass M. One also finds that 

2m + mn + An — m? 1 
%)(4)(r) = - 2 - (19) 

This result is of great importance, because one can show that if 2m + mn + An — m 2 ^ 
than a sufficient condition for the singularity to be strong in the sense of Tipler is satisfied. 

A strong singularity in the sense of Tipler can be defined as follows ||11|| : Let us consider 
non-spacelike geodesic running into the singularity. Then, if the limit of the greatest lower 
bound of the volume element defined by the spatial metric induced by any three independent 
Jacobi fields vanishes as one approaches the singularity along every such geodesic, then the 
singularity is strong in the sense of Tipler. The physical content of this definition is that the 
volume of every extended physical object is compressed infinitely, such that every extended 
object will inevitably be crushed to zero volume. We note, that this is a 'strong' definition 
for destructive singularities, as even weaker singularities would be strong enough to destroy 



any extended physical object. Tipler [JTTJ] gives an example for such a singularity, where the 
spatial metric induced by the Jacobi fields behaves near the singularity like diag(r, r, t~ 2 ). 
The volume element defined by the spatial metric is unity, but any extended object would be 
infinitely stretched in one direction, and infinitely squeezed in two directions, in an infinite 
spaghettification. 



A theorem by Clarke and Krolak |12| relates Ru)U) to the strength of the singularity: 



A sufficient condition for a singularity to be strong in the Tipler sense is that the integral 
Jo* dr' Jq <ir // i?(4)(4)(r") diverges at r*. For example, the FRW cosmological model has a 
Tipler strong singularity, as for this case this integral diverges logarithmically. Note that 
this is a sufficient condition, but not a necessary condition: the Schwarzschild and Kasner 
singularities are Tipler strong, but the integral vanishes identically by merit of their vacuous 
spacetimes. The important conclusion to be made here is as follows: For any non-vacuum 
homogeneous spacelike singularity near which the metric has the power-law form @ the 
singularity is strong in the Tipler sense by virtue of the theorem by Clarke and Krolak and 
by Eq. (0). 

IV. THE UNCHARGED CASE 

In this section we study the spacelike singularity in the presence of a self-gravitating 
scalar field with spherical symmetry and with vanishing electric charge (q = 0). This is 
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the non-linear generalization in the homogeneous case of the linear analysis of Doroshkevich 
and Novikov 0. In the next section we shall also study the charged case. We assume 
a series expansion for the metric functions of the general form f^ n '(r) = J27=i fi( r ) an d 
/i<n)(r) _ Y^ =1 hi(r), where at r = f( n > = h^ 1 ' = 0. Here, n denotes the expansion order 
of the series, which, as will be shown below, are assumed to have a finite radius of absolute 
convergence. That is, j( n=00 ) = f ^ ]i( n =°°) = fa f or some finite interval < r < r$. We 
next assume a power-law series behavior of the metric functions. That is, we assume that 
hi = hio r ni and f\ = fa r mi . Note that for the time being we do not restrict the parameter 
range of rrii. 



A. The leading order approximation 



To the leading order all nonlinear terms in the rhs of Eqs. (|7|) and are negligible, and 
therefore Eqs. (0) and (§) assume the form 

d 2 f 2 

/^-A + ?4 = (20) 



h[r + + — §i = (21) 

Substitution of the Ansatz h\ = B r@, fi = —A r a , with A, B > (as we are looking for 
a spacelike singularity) into these equations yields a = (3 + 2 and A = (/3 + l)B/d 2 . From 
the positivity of A and B it then follows that (3 > — 1. In Schwarzschild, /3 = — 1, and 
the exponents are uniquely determined. Here, however, the scalar field endows the field 
equations with a freedom in the value of (3. For small amplitudes of the scalar field on the 
initial hypersurface we would expect the deviation of (3 from —1 to be small, such that (3 
would still be negative. However, there would also be situations with positive values for (3, 
and even (3 = ClS db special case. We thus find that 

h w =5/ (22) 



f il] = -((3 + l)j 2 r^ 2 . (23) 

Recall that d is a constant of integration for the Klein-Gordon equation @. In Schwarzschild 
d = 0, and therefore one needs to take 

lim = -j- lim B = 2M, 

d ^0 d 2 AM 2 d _^ 

(3^-1 (3^-1 

where M is the mass of the black hole, in order to recover in the limit of vanishing scalar 
field the Schwarzschild solution. An important conclusion to be drawn from this solution is 
that r cx (t>* — v ) 1 ^ 2 . This can be shown as follows: For a general homogeneous metric (0) 
we define a 'tortoise' coordinate r* by g rtrt = —gu- The metric then takes the form 
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ds 2 = g r , r , ( dr 2 - dt 2 ) + r 2 dtf, (24) 
where g TtTt = -g tt = g„ {dr/ dr*) 2 . 



Consequently, dr/ dr* = — J—gtt/grr = — (1/r) w ci 2 /(/3 + 1). Integration yields r* 



— |y(/3 + l)/d 2 r 2 + const. Defining now future-directed null coordinates by t = \{v — u) 
and r* = |(t> + it), we find t> = 2r^ + const, or 

i 

/ iO \ — 

(v*-v) 1/2 . (25) 



/3 + 1 



The power | we find is independent of the value of (3, and is a direct consequence of the 
relation a = (3 + 2. We note that for Schwarzschild a = 1 and (3 = — 1, which clearly satisfies 
this relation. 

We note, that the square-root relation of Eq. (^) can be deduced directly from the 
full partial differential equations. From our numerical results F(u,v) vanishes faster than 
r{u,v) as functions of v — v* near the singularity. Consequently, from Eq. (6) of Ref. |TIj 



it follows that near the singularity (r 2 ) uv Pd 0. This equation can be readily integrated to 
r 2 (u, v) = Xi{u) + x 2 (v ). As at the singularity r vanishes, one finds that along u p = const a 
series expansion yields 



2/ \ / \ d(r 
r [vj = (v — i>*J 



dv 



+ 0[(v - v*) 



and consequently r oc (v — v*) 1 / 2 . 

We next show that our solution is generic. In our solution there are three arbitrary 
parameters, namely (3, d 2 , and B. Without a scalar field the solution is Schwarzschild, 
where there is just one arbitrary parameter, by virtue of Birkhoff's theorem. Therefore, 
with the scalar field, we would expect one additional arbitrary parameter, i.e., two arbitrary 
parameters. Apparently, we have here three parameters. However, the arbitrariness in the 
fixing of B is just a trivial gauge mode, related to the possibility to make an arbitrary 



transformation t —* t' = f(t) (cf. Ref. [|T^] for the analogue in Schwarzschild). Therefore, 
with a scalar field there are only two non-trivial arbitrary parameters, as should be expected. 
Therefore, our solution has the right number of arbitrary parameters, and in this sense is 
generic. Now, for a system of nonlinear equations the notion of a general solution is not 
unambiguous, and there may be other solutions of non-zero measure in solutions space, but 
our solution above also has a non-zero measure, and is therefore generic. 

We can now substitute our solution in the Klein-Gordon equation (|(J), and obtain for the 
scalar field to the leading order in r 

$M = yjp + l lnr (26) 

which diverges logarithmically like the linear scalar field studied by Doroshkevich and 
Novikov 0. We note that the amplitude of the scalar field depends on j3, which is deter- 
mined by the nonlinear equations, whereas in the linear analysis the amplitude in arbitrary. 
Therefore, even the leading order in the expansion for the scalar field is determined by 
nonlinear effects. 
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We now assess the error involved with the consideration of the leading order only of the 
series expansion. For this aim, we consider Eqs. (0)-©- Substituting our leading order 
expression, we find the expressions for the deviation from zero of the rhs of each equation, 
which we shall refer to below as the error associated with the truncated solution at a certain 
order. We find that Eq. (|9|) is satisfied exactly. To the leading order in r 13 the error in Eq. 
(0) is (B 2 /d 4 ) {13 + l) 2 r 2/3+4 , and the error in Eq. <g) is (B 2 /d 2 ) ((3 + 1) r 2 ^ 2 . 



B. The second order approximation 

We shall now find the second order in the series expansions for the metric functions and 
for the scalar field. To the second order the field equations Eq. (|7]) and Eq. (|8|) reduce to 

M 1 - 2 ^)^ 1 -^)^ (27) 

h' 2 r + h 1 f 1 + h 2 + -p = 0, (28) 



where hi and f\ are already known from Eqs. (|22| ) and (|23|). We again assume a power-law 
behavior for f 2 and h 2 , and obtain 

f 2 - (/3 + l) 2 (3/3 + 4) -B 2 

h ~ {(3 + 2) 2 d* [ J) 

h 2 - PiP + 1} * r 2 ^ 2 (30) 
112 ~ {(3 + 2) 2 d 2 r • m 



We thus find that the second-order expansion for / and h is 

-.2/3+4 

1 lP+ U 2 ((3 + 2) 2 d" 

1,(2) _ R p , jj ^ 2/3+2 

+ (/3 + 2) 2 <P 



(31) 



(32) 



We now evaluate the errors, to leading order in r' 3 , in the field equations (0)-© for the 
second-order expansion. For Eq. © the error is 2(B 3 /d 6 ) [((3 + l) 3 (5/3 + 6)/ ((3 + 2) 2 ] r 3/3+6 , 
the error in Eq. (g) is 4(B 3 /d 4 ) [(/? + l) 3 /(/3 + 2) 2 ] r 3/3+4 , and the error in Eq. (|) is 
[03 + l)(/3 2 + 8(3 + 8)/(/3 + 2) 2 ]/rf 4 r 3 ^ 2 . 

The second order expansion we find from the Klein-Gordon equation (|6|) for the scalar 
field is 



ftW /aTTT i j_ V^+l (5/3 2 + 12/3 + 6) 5 ^ 

$ - v^ + 1 lnr + W+W # ' 

Namely, under nonlinear effects the powers of the second (and higher) order term is different 
from the power found in the linear analysis |?J. 
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C. General expression for the series expansion 



We now turn to the general form of the series expansion. Based on the above expressions, 
we seek series expansions of the forms 

oo 

/ = £ a n r^ n a n = a n (d 2 ,(3,B) (34) 

n=l 



5> 

n=l 



„(/3+2)n-2 



b n = b n (d 2 ,P ,B). 



(35) 



Equations (|34]) and (|35|) will solve the field equations (0H© if the following two conditions 
are satisfied: First, the series converge absolutely [this is to ensure that when Eqs. fl5^]) and 
(j35|) are substituted in Eqs. (0)-© the series multiplication theorem will be applicable], 
and second, that the expansion coefficients a n and b n can be found uniquely for any n. If 
these two conditions are satisfied, than Eqs. (0) and (|3~5p represent a generic solution of 
the field equations (as they contain the right number of arbitrary parameters). The first 
condition is hard to prove, without knowledge of the values of the arbitrary functions d 2 , (3 
and B. We thus assume that the first condition is satisfied for small enough values of r. We 
next prove the satisfaction of the second condition. For n > 2, Eqs. ( |34| ) and ([35]) yield the 
following algebraic equations for the expansion coefficients a n and b n : 



d 2 a n + [(/3 + 2)(n + l) - 1] K 



n-l 
k=0 



n—k 



-ih 



(36) 



where 



{[n(/3 + 2) +13 + 1} b x + 2d 2 a x ) a n + {f3+ l)a x b n = c, 



(37) 



n-l 



n—k—1 



c = - \[( n - k )(P + 2) h + d 2 a k a n _ k + b k ajo»-k-i-i \ 

k=l I 1=0 ) 

n-l 

— b\ Y, a l a n-l-l- 
1=0 

These equation will have a unique solution if the determinant of the homogeneous part 
does not vanish for all n. A straightforward calculation yields for this determinant A n = 
n 2 B([3 + 2) 2 ^ (recall that (3 > —1), which proves that the second condition is satisfied. 
The general series expansion for the scalar field is 



P + l lnr + ^> n 



(P+2)(n-l) 



(38) 



71=2 



where c n are the expansion coefficient which can be found uniquely for each order n from 
the Klein-Gordon equation (|6]). We note that in Schwarzschild (3 — — 1, and indeed for this 
choice of (3 the scalar field vanishes identically. 
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V. THE CHARGED CASE 



In this section we shall assume that q ^ 0, and show that there exists a generic solution 
for a spacelike singularity, of similar properties with the singularity one finds numerically 
@||. We shall show that this singularity has many similarities to the singularity in the 
uncharged case we studied in the previous section, but also some important differences. We 
again make the Ansatz that the metric functions are functions of r only, and that r is a 
timelike coordinate. We thus neglect all derivatives with respect to spacelike coordinates. 
We first consider the leading order expression, then the second order correction, and finally 
the general expression for the series expansion. 



A. The leading order approximation 

The source term for the Einstein-Maxwell-Klein-Gordon equations Eqs. contains 
two contributions: the contribution of the electric field and the contribution of the scalar 
field. In order to find the leading order of the series expansion for the metric functions, one 
can consider the following three possibilities for the relative contributions of the two sources 
near the spacelike singularity: 

(a) The contribution of the scalar field near r = is dominant, and the electric field's con- 
tribution is negligible. 

(b) The contribution of the electric field near r = is dominant, and the scalar field's con- 
tribution is negligible. 

(c) The contributions of the scalar field and of the electric field near r = are comparable. 
It turns out that there is no consistent solution of the field equations with possibility (c). 
Let us now consider possibility (b): As near r = the scalar field's contribution to the 
field equations is negligible compared with the electric field's contribution, the leading order 
expression is the same as the leading order expression without a scalar field at all. How- 
ever, from the generalized Birkhoff theorem, the solution then is nothing but the Reissner- 
Nordstrom solution, written as a series expansion near r = 0. But the r = singularity in 
Reissner-Nordstrom is timelike rather than spacelike, which is the type of singularity we are 
interested at here. Therefore, we do not expect possibility (b) to realize in our case. We note 
in passing that this case might be of some relevance for consideration of the self-gravitating 
scalar field solution for an hypothetical extension of the spacetime manifold beyond the 
weakly singular Cauchy horizon. 

We study, then, possibility (a). Namely, we assume that near the singularity the con- 
tribution of the electric field to the energy-momentum tensor is negligible compared with 
the contribution of the scalar field. With this assumption, the equations which govern 
the first order expressions for the metric functions and for the scalar field are the same as 
their uncharged counterparts Eq. Q2(]| ) and Eq. (|2T|). Consequently, to the leading order, 
the spacelike singularities with an electric field or without an electric field have the same 
functional form [cf. Eqs. and (P3"D1: 

=5/ (39) 
/ (1) = -(/5 + l)|r^ 2 . (40) 
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There is, however, a difference in the range of the parameters between the two cases. Whereas 
in the uncharged case (3 was free to assume negative, positive or zero values (subject to the 
restriction (3 > —1), we here find that in order to have a solution consistent with our con- 
clusions based on numerical simulations f3 needs to be positive. Without a scalar field, the 
solution is Reissner-Nordstrom, where the singularity is timelike, and (3 = —2. Because 
there is no scalar field, the electric contribution to the energy-momentum tensor cannot be 
neglected, and the first order equations are dominated by the electric contribution. There- 
fore, there is no smooth analytic limit for vanishing scalar field, and there is a discontinuity 
in the values that (3 can take. This can be understood in terms of the causal structure of 
the singularity, which with a scalar field is spacelike and without a scalar field is timelike. 
With Eqs. (|3~9D and (|4Cf) , the Klein-Gordon equation @ yields for the scalar field again 

$W = y/p + 1 lnr (41) 

which is the same as $W for the uncharged case [Eq. (plf)]. Note that there is no value of (3 
which nullifies the scalar field in this with vanishing scalar field there is no spacelike 

singularity like the singularity we are studying. 

The estimate for errors associated with the leading order terms in the charged case are, 
to the leading order in r p : For Eq. (0) the error is — (B 2 /d 4 ) ((3+ 1) 2 q 2 r w+2 + {B 2 /d 4 ) ((3 + 
l) 2 r 2/3+4 , for Eq. (§) the error is -(B 2 /d 2 ) (f3 + l) q 2 r 2/3 + (B 2 /d 2 ) ((3+1) r 2/3+2 , and for Eq. 
(§) the error is 2(B 2 /d 2 ) ((3 + 1) q 2 r 2/3 ~ 2 . Note the differences between these error estimates 
and their counterparts in the uncharged case. For each equation, the error estimate here 
has two terms, one of them does not depend on the charge, and the other does depend 
on it. The term which does not depend on the charge is the same as the error term in 
the uncharged case. However, the charge- dependent term has a smaller power index, and 
is therefore dominant. [In the error for Eq. @ the term which does not depend on the 
charge vanishes.] The difference in the power indices between the charge-dependent and 
the charge-independent terms is 2. We shall show below that this is an indication for a 
more complicated series form for the metric functions than the form we encountered in the 
uncharged case. 



B. The second order approximation 

The electric field does not contribute then to the leading order terms in the series ex- 
pansion for the metric functions or the scalar field. However, we find that it does contribute 
to the second order terms. The equations which govern the second order terms are: 

h 2 r + h 2 + —3- + q —y = 0, (43) 

where again fx and f 2 are known from the first order terms ( [39|) and (0), and where we 
again assumed a power-law behavior. These equations are different from their uncharged 



14 



counterparts [Eqs. (|27|) and (^)] not only in the appearance of charge-dependent terms, 
but also in the absence of terms which appeared in the uncharged case. Substitution of the 
power-law behavior Ansatz yields 

r = (g+lgjg+2) & , 2/3+2 

J2 /5 2 d 4 q 1 J 



09 + 1)09 + 2) B 2 ^ 2/3 



& 1 



IP' T 



(45) 



such that 



(2) 



(8 + 1)* rf»* + (^+l) 2 (3^ + 2) &_ , rW+2 



(46) 



hV=Br*-V + W + V*f 1 *' (47) 
and from the Klein-Gordon equation @ 

*<» = #TT in, - ^y^' £ 9 * * (48) 

The estimate for the leading terms in the errors associated with the second order terms 
are: For Eq. (0) the error is ((3+1) 2 (B 2 /d 4 ) r 2/3+4 + 2[(/3+l) 3 (5/?+4)//5 2 ] (B 3 /d 6 ) q 2 r 3(3+2 + 
0(r 3/3+4 ), for Eq. (@) the error is (/3 + 1) (5 2 /rf 2 ) r 2/3+2 + 4[(/3 + l 3 )//? 2 ] (5 3 /d 4 ) g 4 r 3/3 + 
0(r 3/3+4 ), and for Eq. (g) the error is + 1) 2 (/3 2 — 6/3 — 8)//5 2 ] (5 3 /d 4 ) g 4 r 3m - 2 + 0(r 4/3 - 2 ). 

Comparing these expressions for the errors associated with the second order expansions 
and the errors associated with the first-order expansions reveals that the effect of the second- 
order terms was to provide terms which balance the leading-order errors for the first-order 
expansions, and leave just the subsequent terms of the first-order expansions as the leading 
terms for the errors of the second-order expansions. Namely, the contribution to the solution 
of each order balances the leading-order terms in the errors associated with the lower-order 
solution, and thus reduces the errors. 



C. General expression for the series expansion 

One would be tempted to consider a series form similar to the series form we have for 
the uncharged case (|34j) and fl35|). However, when one considers the third order terms for 
the metric functions fs and /13 in the charged case one finds that this is impossible, unless 
the exact value of (3 is known. A close inspection of our previous results and of the third 
order approximation shows that the increments in the power indices in the charged case are 
not equal. Namely, the increment changes from /3 to (3 + 2. Therefore, it is more natural to 
seek a general solution in terms of a double series expansion than a single series. Therefore, 
the solution will be of the form 
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/ = E E amn ^ (m - 1)+(/3+2) " a mn = a mn (d 2 , g 2 , 5, /?) (49) 



m=l n=l 



oo oo 



ft=EE & ™ r /3m+(/3+2)(n-l) ^ = ?2> s> ^ (5Q) 



m=l n=l 



Again, in order that Eqs. (|49"1) and (p0|) would be a solution of the field equations 
the series should converge absolutely for some finite convergence interval, and the series 
expansion coefficients a mn and b mn should be found uniquely for each m,n. In fact, the 
terms we previously found explicitly are nothing but the an, 031 and fen, 621 terms in this 
double series expansion. It is still an open question whether the satisfaction of these two 
conditions can be proved rigorously A similar double series expansion is also obtained for 
the scalar field, with the leading term diverging logarithmically. 

VI. SINGULARITY STRENGTH 

From the expression for the metric functions we found for the charged case we can 
calculate the leading term for the Kretschmann scalar and for Ru)U) in terms of the proper 
time of an observer who follows at — const radial trajectory. We find that 

Ra ^ R ~ 64 (/? + 4)4 J^rf (51) 
%)(4)W=8t^3' (52) 



(/3 + 4) 2 r 



From Eq. (^) and from the theorem by Clarke and Krolak |12| it then follows that the 
spacelike singularity in spherical charged black holes with a self-gravitating scalar field is 
strong in the Tipler sense. Recall that the singularity we found is a generic one, in the sense 
that it relies on the correct number of arbitrary parameters. One indeed finds that these 
are the expressions which one obtains from the general expressions for power-law metrics 
(|18|) and (|l9l) for this solution with m = (3 and n = (3 + 2. In addition, as in the charged 
case (3 > 0, it turns out that there is no special case which nullifies either of the curvature 
invariant Ra^sR ^ 6 or the second integral over i?( 4 )( 4 )(r). Namely, for any choice of (3 
the spacelike singularity is strong in the Tipler sense. (Note that because of the inevitable 
presence of the scalar field for the singularity to be spacelike with electric charge, there is 
no vacuum solution, and therefore the Ricci tensor does not vanish.) 



VII. CONCLUDING REMARKS 

We found a generic solution for a spacelike singularity for the spherically symmetric 
Einstein-Maxwell-scalar field equations, which was previously found numerically in scalar- 
field collapse simulations. The generic singularity we found has the same properties as the 
singularity which arises in the numerical simulations. 
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The singularity we found was obtained under the assumption of quasi-homogeneity, 
namely, that the metric functions do not depend on derivatives with respect to the spa- 
tial coordinates (in our case, because of the spherical symmetry, this means that there is no 
dependence on t). Namely, we sought a velocity-dominated singularity. The next obvious 
step in the analysis of the singularity would be to relax the homogeneity assumption. This 
could perhaps be done by allowing the arbitrary parameters in the solution to be functions 
of t. Namely, all we need to do is to make the transformation /3 — > [3(t), B — > B(t), and 
d 2 — > d 2 {t) for the leading order terms of the solution. The higher-order terms will have 
to be corrected for the terms with dependence on derivatives with respect to t. (Because 
of the spherical symmetry there would still not be derivatives with respect to the angular 
coordinates 9 and 0.) A preliminary check shows that these leading terms indeed satisfy the 
full spherically-symmetric Einstein-Maxwell-scalar field equations, with error smaller than 
the leading terms. However, it is still needed to perform a fuller analysis of the higher-order 
contributions. Then, the solution would be inhomogeneous, and depend on three arbitrary 
functions, and in this sense would be a generic solution. The homogeneous solution is then 
nothing but the pointwise behavior of the inhomogeneous solution. However, it is still to be 
shown that this inhomogeneous solution really solves the inhomogeneous Einstein-Maxwell- 
scalar field equations. 

Special attention is required in order to analyze to behavior of the fields near the space- 
time event corresponding to the change of the causal structure of the singularity, namely 
to the event where the null singularity becomes spacelike. There are two possibilities to 
approach this spacetime event numerically, using a double-null code such as ours. First, 
one could approach this event along an outgoing null rays. Then one would need to reduce 
the value of the outgoing null coordinate u, and compare the values of the fields approach- 
ing the singularity along various values of u. Second, one could also study the singularity 
along ingoing null rays. Late-time ingoing null rays actually probe the null mass-inflation 
singularity. However, one can also focus attention to the deep portions of these rays, and 
thus study the nearly-complete focusing domain, where the value of r is very small. We 
find numerically that the area coordinate r continues to decrease monotonically even at 
the nearly-complete focusing domain. This second approach would be more natural with a 
sleight variation of the numerical code: The integration is normally carried out in our code 
along ingoing rays, namely, we intergate from the initial data on an ingoing ray and the 
data on the first grid-point on the next ingoing ray, the fields throughout this second ray, 
and so on (for more details see Ref. [TIJ). However, for this second approach to the special 
event where the causal structure of the singularity changes, one could benefit from changing 
the direction of integration, and integrate along outgoing null rays instead of ingoing rays. 
Then, for each value of advanced time v for a certain ingoing null ray one would eventually 
approach the spacelike singularity for large u, and the later v , the closer one would be the 
the spacetime event under consideration. 

We again stress that it is still an open question whether in a more realistic model our 
results would be preserved. Namely, a scalar field was introduced as a toy model for a 
physical field because it has a radiative mode in spherical symmetry. However, a scalar field 
is not a realistic physical field, and it is possible that more realistic fields will not create a 
spacelike singularity. It has been shown before that scalar fields create unique phenomena 
(see, e.g., Ref. 0), and therefore it is not impossible that the generic spacelike singularity 
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we found is just an artifact of an over-simplified toy model. In addition, a realistic black 
hole will not be spherically symmetric, and whether a spacelike singularity would be created 
inside spinning black holes is a question still awaiting investigations. 
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